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Abstract 

An outstanding problem in statistical mechanics is the order parameter 
of the chiral Potts model. An elegant conjecture for this was made in 1983. 
It has since been successfully tested against series expansions, but there is 
as yet no proof of the conjecture. Here we show that if one makes a certain 
analyticity assumption similar to that used to derive the free energy, then 
one can indeed verify the conjecture. The method is based on the "broken 
rapidity line" approach pioneered by Jimbo, Miwa and Nakayashiki. 
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1 Introduction 

The chiral Potts model was originally formulated as an A-state one-dimensional 
quantum hamiltonian ^ [2] , and then as a two-dimensional classical lattice 
model in statistical mechanics. |IJ [3] It satisfies the star-triangle relations. 
The free energy was first obtained for the infinite lattice using the invariance 
properties of the free energy and its derivatives. [H] Then in 1990 the functional 
transfer matrix relations of Bazhanov and Stroganov [Jj were used to calculate 
the free energy more explicitly as a double integral. [51 l9l ITU] 

The chiral Potts model is a system of spins living on the sites of a planar 
lattice, usually (but not necessarily ^JJ) taken to be the square lattice C. 
Each spin takes one of N possible states, labelled 0, . . . , N — 1 and interacts 
with its neighbours. If the spin at site i is (jj, then the interaction between 
spins on adjacent sites i and j depends on crj, aj only via their difference 
<jj — (jj (mod N). 

We shall define the interactions in the next section. Here we merely note 
that they depend on a temperature- like parameter k' , which is small at low 
temperatures and large at high temperatures. For k' < 1, a related parameter 
is 

k = y/(l- k' 2 ) . (1) 
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Let oj = e 2m / N and let a be a spin deep inside the lattice. Define 

Mr = (to ra ) (2) 

as the average value of uj ra , for r = 0, . . . , N. Suppose one fixes the boundary 
spins to be zero and allows the lattice to become infinitely large, a remaining 
near the centre. Then for sufficiently high temperatures (weak enough inter- 
actions) the boundary conditions will become irrelevant, all values of a will 
be equally likely, and M r will be exactly zero for r = 1, . . . ,N — 1. 

However, the system displays ferromagnetic order. There is a critical tem- 
perature T c below which the boundary conditions remain relevant even for 
an infinitely large lattice. Then Mi, . . . ,Mn-i are non-zero. They can be 
thought of as spontaneous magnetizations or order parameters. We expect 
M r to vanish as T — > T c , being then proportional to 

(1-T/T c f , 

the index [5 being known as a critical exponent (dependent on r). 

In fact the critical point is when k' = 1 and k = 0, so the ferromagnetic 
region is when < k, k' < 1. For N = 2 the chiral Potts model reduces to the 
Ising model. From the exact results of Yang |13| and Onsager|14j we know 
that then Mi = fc 1 / 4 . Hence Mi = (1 — A;' 2 ) 1 / 8 and the critical exponent (3 
is 1/8. 

In 1983 Howes, Kadanoff and den Nijs P considered the case N = 3 and 
evaluated Mi,M^ to order A;' 13 in a series expansion in powers of k 1 . They 
found their series fitted the formula Mi = M2 = k 2 ^ 9 , giving j3 = 1/9. Later, 
Henkel and Lacki ^S] expanded ^2M r for general N to order k' G . 

In 1989 Albertini et al jTHJ expanded the general N case to order k' 5 , and 
found that all the results were consistent with the remarkably simple and 
elegant conjecture 

Mr = k r( ^ N - r ^ N2 , 0<r<N (3) 

(see also Ref. 1 1 7 j ) . giving (3 = r (N - r)/2N 2 . Baxter [TT] used finite corner 
transfer matrices to expand M r for N = 3 to order &/ 14 , and again found the 
results were consistent with the conjecture (|3*j). 

So for the last fifteen years or more there have been conjectures for the 
order parameters of the chiral Potts model. However, the author knows of no 
derivation of any of these for N > 2. One can contrast this with the situation 
for the Ising model, where it was five years from the time Onsager calculated 
the free energy to when at a conference in Florence in 1949 he announced 
the formula for Mi, and three years from then until Yang published a proof. 
(The actual calculation took Yang six months of work off and on |18t p. 12].) 

The difficulty has been that unlike most other solvable models, the chiral 
Potts model does not have the "rapidity difference property" . For this reason 
the infinite-size "corner transfer method" for exactly calculating the order 
parameters fails. An alternative method was developed in 1993 by Jimbo et 
al ^Hl- It relies on breaking one of the rapidity lines of the lattice, giving one 
half- line a different rapidity value from the other. The author applied this 
method to the chiral Potts in 1998 [20] and wrote down functional relations 
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for the resulting generalized order parameters G pq {r). They are functions of 
the rapidities p, q of the two half-lines. 

Again there was a difficulty. The relations by themselves do not completely 
define the G pq {r). One also needs information as to the analyticity properties 
of the functions. 1 In this respect the functional relations are similar to the 
"inversion relations" for the infinite-lattice free energy |211 I22| . 

The calculation of the free energy of the chiral Potts model [9|ll0|l22] pro- 
ceeds in two stages. First one considers a related il T2(t p )" model. This is in- 
timately connected with the superintegrable case of the chiral Potts model. [23] 
It is much simpler than the chiral Potts model in that its Boltzmann weights 
depend on the horizontal rapidity p only via a single parameter t p , and are 
linear in t„. Its row-to-row transfer matrix is the product of two chiral Potts 
transfer matrices, one with horizontal rapidity p, the other with a related 
rapidity p(0, 2) defined in eqn. (|12j) of section 3. 

For a finite lattice, the partition function Z of the T2(t p ) model is therefore 
a polynomial in t p . The free energy is the logarithm of Z^ M , where M is the 
number of sites of the lattice, evaluated in the thermodynamic limit when 
the lattice becomes infinitely big. This limiting function of course may have 
singularities in the complex t p plane. A priori, one might expect it to have N 
branch cuts, each running though one of the N roots of unity. However, one 
can argue that in fact it only has one such cut. As a result the free energy 
(i.e. the maximum eigenvalue of the transfer matrix) can be calculated by a 
Wiener-Hopf factorization. 

The second stage is to factor this free energy to obtain that of the chiral 
Potts model. 

It turns out that the first stage of this process can be applied to the 
generalized order parameter function G pg (r), provided we take q, p to be 
related by eqn. (|B2*|) of section 6. 

We present the working in the following sections. We do not need to 
calculate G pq (r) for general p, q and we do not do so. 

We define the model in section 2, and the function G pq (r) in section 4. We 
also present the functional relations satisfied by G pq (r), but in fact we hardly 
use them. It is the analyticity properties that hold the key to calculating 
the G pq (r), and most of the discussion in sections 3 to 6 is concerned with 
presenting evidence for our assumptions regarding these properties. 

In particular, we show that when p, q satisfy (jri2j) . the G pq (r) can be ex- 
pressed in terms of T2{t p ) Boltzmann weights U(a, b, c, d) that are linear in t p . 
We argue that G pq (r) is therefore like the free energy of the T2(t p ) model, in 
that it has at most one branch cut in the i p -plane, rather than the N cuts 
that one might expect. 

We give our precise assumptions in section 7. Then in section 8 we use 
them to obtain G pq {r) by a Wiener-Hopf factorization in very much the same 
way as one calculates the T2(t p ) free energy. The desired formula © follows 

1 The example I like to quote is the relation f(z + 1) = f(z). By itself this merely says that 
/(z) is periodic. However, if one can also show that f(z) is analytic and bounded in the strip 
< 'St(z) < 1, then it follows that f(z) is analytic everywhere and bounded, so by Liouville's 
theorem it is a constant. 
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immediately. We also present an alternative method, looking at the product of 
N such functions G pq (r), that avoids the need for a Wiener-Hopf factorization. 

In section 9 we briefly discuss some other special cases (analogues of the 
Tj(t p ) models for j = 3, . . . , N) that may be tractable, and make a conjecture 
for the form of G pq (r) for some of these cases. 



2 Chiral Potts model 

We use the notation of [SlIHlEnj- Let k,k' be two real variables in the range 
(0, 1), satisfying Consider four parameters a p ,b p ,c p , d p satisfying the ho- 
mogeneous relations 

dp ~t~ k bp = kdp , 

k'a^ + bp* = kc^ . (4) 

Let p denote the set {a p , b p , c p , d p }, or rather their ratios, since all the equa- 
tions we shall write involve a p ,b p ,c p ,d p only via their ratios. Similarly q 
denotes the set {a q ,b q ,c q ,d q } satisfying the relations ® (with p replaced by 
q). We call p and q "rapidity" variables. 

Define Boltzmann weight functions W pq (n), W pq (n) by 

= n f;- a ^i , 

b p d q - c p a q uj3 



where 



J = 1 C p b q - b p CqUj3 



Then the conditions Q ensure that 

W pq (n + N) = W pq (n) , W pq (n + N) = W pq (n) , 
so the functions are periodic of period N. Note that 

Wpq(0) = W pq {0) = 1 . 

Define related parameters 

x p — dpi dp , y p — bp/ Cp , tp — XpUp , fip — dp/ Cp . (6) 
They satisfy 

x£ + = k(l + x?y?) , kx» = l- k'/fj,p , ky? = 1 - , (7) 

kX = (1-A/<)(1 -*'/<) . 
We can also write © as 



W pq (n) = II , (8a) 

j—l Up ~ "J-* JL q 
n j 

W pq {n) = M» II ^ 'Y" ■ (8b) 
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Fi gure 1: The square lattice C of solid lines and circles, showing the central spin a and the dotted 
directed rapidity lines . 



Now consider the square lattice C, drawn diagonally as in Figure ^ with 
a total of M sites. On each site i place a spin G{, which can take any one of 
the N values 0, 1, . . . , N - 1. With each SW - NE edge (with % below 

j) associate an edge weight W vq {pi — Oj). Similarly, with each SW - NE edge 
(j; k) (j below k), associate an edge weight W pq (aj — <7fc). Then the partition 
function is 

Z = E II Wnfa - aj) n Wpgfo - <r k ) , (9) 

the products being over all edges of each type, and the sum over all N M values 
of the M spins. We expect the partition function per site 

to tend to a unique limit as the lattice becomes large in both directions. 

Let a be a spin on a site near the centre of the lattice, as in Figure^ and 
let /(a) be any function thereof. Then the thermodynamic average of f(a) is 



(/(a)) = r 1 ^ /(a)n^-^)II^-^) 



(10) 



we expect this to also tend to a limit as the lattice becomes large. 

We have also shown in Figure^the vertical and horizontal "rapidity lines" . 
Each edge of C passes through the intersection of two rapidity lines. 

We could allow each vertical (horizontal) rapidity line a to have a different 
rapidity p a (q a ). If an edge of C lies on lines with rapidities p a , qp, then the 
Boltzmann weight function of that edge is to be taken as W pq (n) or W pq (n), 
with p = p a and q = qp. 

The weight functions W pq (n), W pq (n) satisfy the star- triangle relation. 
For this reason we are free to move the rapidity lines around in the plane, 
in particular to interchange two vertical or two horizontal rapidity lines. [25] 
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So long as no rapidity line crosses the site with spin a while making such 
rearrangements, the average (/(a)) is unchanged by the rearrangement. 2 

All of the rapidity lines shown in Figure ^ are "full" , in the sense that 
they extend without break from one boundary to another. We can move any 
such line away from the central site to infinity, where we do not expect it to 
contribute to (/(a)). Hence (/(a)) must be independent of all the full-line p 
and q rapidities. 

It can still depend on k or k', which play the role of universal constants, 
the same for all sites of C As we mentioned in the introduction, it has been 
conjectured that 

Mr = (w ra ) = k ri - N - r ^ N2 , 0<r<iV . (11) 

The aim of this paper is to verify this conjecture, subject to a plausible ana- 
lyticity assumption. 

3 7*2 model 

One of the difficulties of the chiral Potts model is the multi-valued nature of 
the relations between x p ,y p ,fi p . Every such relation involves taking an iVth 
root. For N = 2 the model reduces to the Ising model. In this case there is a 
simple uniformizing substitution whereby all variables can be written as single- 
valued Jacobi elliptic functions of another variable u r . \22l App. B] For higher 
values of N no such substitution is known. It is therefore significant that 
the chiral Potts model can be related to "superintegrable" or u Tj(t p )" models, 
where the dependence of Boltzmann weights on the horizontal rapidity p is 
simple: they are explicit polynomials in the single variable t p . (They still 
involve all the vertical rapidity variables x v ,y v ,fi v .) 

The key equations are given in [H]. Let us change notation and use the 
symbols v or v' for the vertical rapidities, p, q for the rapidities of the lower 
and upper rows. Thus we replace the p,p', q, r of [H] by v, v',p, q. 

For generality, we allow v and v' to de different in this section. In later 
sections we take v' = v and all vertical rapidities to be the same, p and q will 
play the role of variables in the functional relations we discuss, while we shall 
regard v as a constant. 

Consider two rows of edges of C Let the horizontal rapidity of the lower 
row be p, and of the upper row 

q = P(k,£) , (12) 

by which we mean 

x q = u) k y p , y q = uj e x p , fi q = 1//j, p , (13) 
k, I being integers. We can impose the restriction 

1 < k + t < N . (14) 

2 Subject to boundary conditions: here we are primarily interested in the infinite lattice, where 
we expect the boundary conditions to have no effect on the rearrangements we consider. 



6 



If we sum over an intermediate spin between these rows, we construct a 
combined weight function 

N-l 

U(a,b,c,d) = Y. W vp( a -9)W vq {g-d)W v , p {b-g)W v , q {g-c) , (15) 

9=0 

depicted in Figure [21 This can in turn be written as 

N-l 

U(a,b,c,d) = N ^ V vpg (a, d; n) V v > qp (-c, -b; n) , (16) 

n=0 

where 

Af-l 

V vpq (a, d; n) = N' 1 ]T W vp (a - g)W vg (g - d) . (17) 

9=0 




Figure 2: The combined weight U(a, b, c, d) of four edges, summed over the middle spin g. 

Let C,ki be the set of integers {—k, 1 — k, ...,£— 1}, with k + I elements. 
We say that n e Cki if n is equal, modulo N, to one of the elements of Q^t- 
Then in (3.16), (3.17) of |Hj it is shown that 

V vpq (a, d;n) = if a - d e (k£ and n i Ck£ , 
V vqp (-c,-b;n) = if n e C,U and b-c<jt( ke , (18) 
U(a,b,c,d) = if a — d e Cki and b — c £ Cki ■ 

Consider a row of adjacent spins ai, . . . , cr s , below a similar row a[, . . . ,a' s , 
as in Figure 01 with intervening face weights (|15j) . For the moment ignore the 
spin g and the dotted lines. The combined weight of the rows is 

s-l 

I(a u ...,a' s ) = n^i^i+lX+i,^) (19) 

i=l 
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Figure 3: Two rows of spins with face weight functions U(ai,<Ji + i,a' i+1 ,a' i ) in between, and a 
spin g to the left linked by edge weights to o\ , a[ . 



and it follows that if o\ — a[ e Chj then the weight is zero unless a% — e Cfc£ 
for i = 1, . . . , s. In this case we only need U(a, b, c, d) for a — d and b — c in 
Cm. From (3.21), (3.39) of 8,, if we define 



then 



j = k + t , a = mod (a — d + k, N) 
m = mod (n + k,N) , /5 = mod (b — c + k, N) , 

V vpq (a,d;n) = fi« uj nd ~ mk (b p /d p ) a y- m f vp (j, a, m) 



(20) 
(21) 



V v/qp (-c,-b;n) = n k J p h%LJ kc - mb (d p /b p ) 



T]p,j,m 



provided a — d,n,b — c are all in which means that 

< a, m, (3 < j < N . 



(22) 



The expressions (|21[) become less frightening if one groups them into fac- 
tors of various types: 

(i) Factors independent of a,b,c,d,n,a, /3,m. 

These are the factors Q kl p , Qyl p , bf$L. They are defined in (3.24) and (3.35) 
of [H]. In this paper we do not calculate full partition functions, but rather 
ratios of partition functions such as (|34|) and (|47|). These factors cancel out 
of such ratios, so we can ignore them. Hereinafter we shall simply take Q kl p = 

T) kl - h ij) 



1 in (|2]J. 
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(ii) Factors that are powers of {b d /d p ). 

These are the factors (b p /d p ) a , (d p /b p )P. In Figure 01 each internal vertical 
edge (aj,(Tj) acquires a weight (d p /b p )P from the face on its left, and a factor 
(d p /b p )~P from the face on its right, where (5 = mod (aj — a'j + k,N). The 
contributions to the internal edges therefore cancel. 

(iii) Factors y~ m , y™. 

These cancel from the product on the rhs of (|15|) . 

(iv) The oo, 7], f vp factors. 

These remaining factors depend on the rapidity p only via the functions 
Vp,j,a, fvp(j,Oi,m). These are defined in (3.26), (3.37) and (3.38) of jS] (where 
our f vp becomes F pq ). They have a vital property: they are Laurent polynomi- 
als in t p , with no other dependence on p. More strongly, rjpj^a is proportional 
to t p , and f vp (j, a, m) is a polynomial in t p of degree not greater than m and 
j — a — 1, with minimum power not less than zero and m — a. 

It follows from d ) that (dp/bp) 01 '? U(a, b, c, d) is a polynomial in t p of 
degree j — 1. Hence from (j!9|) 

I(ai,...,a' s ) = {bp/dp) x ~»l{o- u ...o-' s ) , (23) 

where T(o\, . . . a' s ) is a polynomial in t p of degree (s — — 1) and A = 
mod (cji — cr[ + k, N), v = mod (a s — a' s + k, N). 
A related quantity that we shall need is 

J r (ai,...a' s ) = V vlqp {-a' 1 ,-ax;-r)I{a l ,...a' s ) , (24) 

where r is some integer. This corresponds to including the spin g in Figure 
13 together with weight functions between g and a\, a[, multiplying by u/ 9 , 
and summing over g. The type (ii) factors now cancel from the left-hand edge 
(ai,a[), but a type (iii) factor arises from the extra V v / qp term. Using 
Vp = tp/xp, it follows that 

xf r U r (a 1: ...a' s ) = {dp/bpYtp- 1 J r {a l ,...a' s ) , (25) 

where m(r) = mod (k— r, N), v = mod (a s — a' s + k, N), and J(t p ) is a Laurent 
polynomial in t p . 

The case j = 2 

One case in which we shall be particularly interested is when 

k = , 1 = 1 , j = 2 . (26) 

For this case the functions 77, f vp (j, a, m) are given in (3.48) of [2|. From (|16|) . 

N ( b \ a+d ~ b ~ c 

U(a,b,c,d) = ^ U(a,b,c,d) (27) 

VvVv' \dpj 
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a — d 


b-c 


U (a, b, c, d) 








y v y v , - u d - b+1 t p 





1 


-ufi v d p (y v - u; d ~ 6+1 av) 


1 





Hv{Vv' - uj d ~ b+1 x v ) 


1 


1 


—UJfl v fl v /(t p — u) d ~ b+1 x v x v i) 



Table 1: The face weights U(a,b,c,d) of the T%(t p ) model. 

where U(a, b, c, d) is given in Tabled 

In particular, from l|19jl and Q23|) it follows that 

1(0, . . . , 0) = (y v y v , - uj^+HpY- 1 , (28) 

ignoring factors independent of p. 

The case j = N 

The other case we shall need is when 

k = -l , £ = N + 1 , j = N . (29) 

Then (^i is the full set of N integers 1, . . . , N and V vpq (a, d; n), V v>qtP (— c, —b; n) 
are always given by (|21[l. with a, (3, m in the range [0, N — 1]. 

In this case we shall be interested in the modified product (|24|) . when 
a s = a' s = 0. Then J r (a\, . . . , a' s ) is given by (|25|). with i/ = JV — 1. The factor 
(dp/bp) u is now independent of r and the spins <ti, . . . ,cr' s _ 1 . Like the type 
(i) factors, it cancels out of the ratios of interest (|34|) and (|47|) . The function 
Jripi-, ■ ■ ■ cr's) in Q25|) is a polynomial in t p of degree (s — 1)(N — 1). 

One particular sub-case that we shall consider is when a\ = ■ ■ ■ = a' s = 0. 
From (3.26), (3.37), (3.38) of |Hj we can verify that 

V P ,N,a = tp . (30) 

f vp (N,N-l,m) = {d v /b v ) N ~ l x™ , (31) 
and hence from (|15|1 and (|21j) that 

AT-l 

[7(0,0,0,0) = (x„x„/ - ^'tp) 5 " 1 , (32) 
i=i 

ignoring factors independent of p. It follows that 

AT-l 

J r (0,...,0) = l[(x v x v , -uhp) 8 - 1 . (33) 
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4 The generalized order parameter 



Now let us replace all the vertical rapidities p in Figure Q by v. We also 
replace all horizontal rapidities q by h, except for the one line immediately 
below the spin a. Following Jimbo, Miwa and Nakavashiki|19j . we break this 
line immediately below the site i containing the spin a, and give the half-line 
to the left the rapidity p, that to the right the rapidity q, as indicated in 
Figure 0] Let F pq {a) be the probability that the spin at site % is in state a, 
i.e. from 

F pq (a) = (S^a) = Z(a)/Z , (34) 

where Z(a) is the sum-over-states with spin cij fixed to be a, divided by the 
unrestricted sum Z = Z(0) + • • • + Z{N — 1). 
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Figure 4: First picture of the function F pq {a). 



Because of the star-triangle relation, F pq (a) is independent of the back- 
ground rapidities v and h, but it does depend on p and q. This is because the 
ends of the half-lines near a cannot be moved away from a. 

In [20] we show that one can obtain functional relations satisfied by F pq (a). 
Let R, S be the automorphisms defined in jS] : 

{a-Rp, bR p , cr p , d,R p } = {b p ,Lua p ,d p ,c p } , (35) 

{asp,bsp,cs p ,ds p } = {uj~ 1/2 c p ,d p ,a p ,uj~ 1/2 b p } , (36) 
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so 

XR P = y P , yn P = ujx p , /j, Rp = l/fi p , (37) 
xs P = ^-/y P , ys P = l/x P , Lis P = u~ 1/2 y p /{x p fi p ) . (38) 
Then from equations (18), (19) of 0, 

W pq (n) = W qtRp (-n) , W pq (n) = W qtRp (n) 

W Rp>Rq (n) = Wpq(-n) , W Rp>Rq (n) =W pq (-n) (39) 
Ws P ,sq(n) = W qp (n) , W Sp ,Sq{n) = Wq P (-n) 
W Sp ,RSq(n) = Wpq(-n) , W Sp ,RSq(n) = Wp q (-n) . (40) 

Symmetries 

From ()39|) . operating by R on all rapidities is equivalent to negating all spins. 
Hence 

F Rp , R q(a) = F p q(-a) . (41) 

There is a reflection symmetry that is not given in .20, • Operate by S 
on all the vertical rapidities, and by RS on all the horizontal rapidities in 
Figure 0] Then from (|40j) this is equivalent to interchanging the functions W, 
W, and negating all spins. Interchanging W with W is in turn equivalent to 
mirror-reflecting the lattice about the central vertical line thorough a, while 
interchanging p with q. Thus 

F R s P ,RSq(a) = Fq P (-a) . (42) 

The dotted rapidity lines in Figures ^ and are directed, bearing arrows 
that give their direction. They form a graph Q (the square lattice) of coor- 
dination number 4. The dual of G is a bi-partite graph, and one of the two 
sub-graphs is the lattice L. 

Although one cannot remove the ends of the half-lines from the spin a, 
one can rotate them subject to the rules given in |20| . notably that one is 
not allowed to introduce any directed circuits into Q. The effect of this is to 
deform Q, but the sites of C continue to live on one of the two sub-lattices 
dual to Q. Every edge of C passes through the intersection of two dotted 
rapidity lines of Q. If the two arrows on Q lie on either side of the edge , 
with rapidities p, q, oriented as for the edge (J, k) of Figure then the weight 
function is W pq (k — j). If both arrows lie on one side of the edge, oriented as 
for the edge then the weight function is W pq (k — j). 

We show in Figure 8 of P0 that if the arrows of all rapidity lines crossing 
a given line of rapidity p point to the left (right), one may reverse the arrow 
on p and replace p by Rp (R~ 1 p). 

The result is that one can perform the following sequence of operations: 

a) rotate the left half-line p clockwise through 90° to a vertical position 
below a, pointing downwards 

b) replace p by R _1 p and reverse the arrow to point upwards, 

c) rotate this half- line R~ 1 p clockwise though another 90° to the horizontal 
position of Figure [3 
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Figure 5: Second picture of the function F pq (a). 



The probability F pq (a) is unchanged by this sequence, so we can alterna- 
tively use Figure [5] to define F pq (a). We can interchange the rapidity half- 
lines q,R~ 1 p therein without changing F pq (a), so F pq (a) = F p / q i(a), where 
q' = R~ l p, R~ l p' = q. Hence 

Fpqifl) = FRq,R-i p (a) , (43) 

which is the relation (15) of [20] . 

Alternatively one can similarly rotate and reverse the left half-line anti- 
clockwise round and above a to lie above q, with p replaced by Rp and again 
with both arrows pointing to the left. This is the configuration of Figure El 
It too can be used for the definition of F pq (a). 

Now if we pass the half-line Rp under q we obtain a different figure, where 
a is connected to the other sites only via a single edge of £, to a site b to 
the right of a. This edge passes though the intersection of the half lines Rp, q 
and has weight W R P , q {b — a). If we remove it, we regain FigureEJ but with q 
interchanged with Rp and a replaced by b. Hence 

N-l 

F Pq {a) = C pq J2 W Rp,'i( b - a ) F R- 1 QM b ) > ( 44 ) 

b=0 

where £ pq is independent of a. With a slight change of £ pq , this is the relation 

(i6) of nag. 
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V V V V V V 



Figure 6: Third picture of the function F pq (a). 

Domains 

If x p ,x q ,y p ,y q ,ujx p all lie on the unit circle in an anti-clockwise ordered se- 
quence, then the functions W P q(n),W pq (n) are real and positive. This is the 
case we have in mind in this paper: it ensures that k and the generalized order 
parameters G pq (r) exist in the thermodynamic limit of a large lattice and are 
continuous and infinitely differentiable. 

For other values of p, q we define the functions to be the analytic continu- 
ations from this physical regime. To fix our ideas, it is helpful to consider the 
low-temperature case, when k' is small. If fi p is of order k', then x p is free to 
take most values in the complex plane, other than those near the N roots of 
unity 1, u>, . . . , w^ -1 , while y p has to be near such a root of unity. If [i p is of 
order k'~ , then it is y p that is free, x p that is constrained. In the Tabled we 
show four such "domains" . They form part of an increasing sequence: if p is 
in domain T> r , then Rp is in domain T> r+ \. 

We also introduce the concept of a "bridge" B r linking adjacent domains 
V r and V r+ \. On such a bridge /j, p is greater than O(k') and less than 0(A;' _1 ), 
and x p , y p are constrained by both V r and V r+ i. For instance, if p is on the 
bridge £>i, then x p ~ 1, y p ~ u. 

If p is in V T and q is on the bridge B r , then W pq (n), W pq (n) are small for 
n / 0. The same is true if p is on B r and q is in T> r+ \. The dominant con- 
tribution to the partition function is then from all spins being zero. One can 
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Domain 


Xp 


y P 






to- 1 


free 


Oik'- 1 ) 


v 


free 


1 


O(k') 




1 


free 


Oik'- 1 ) 


v 2 


free 


CO 


O(k') 



Table 2: Some physical rapidity domains 



develop the usual series expansions for n for any such case: the results are all 
analytic continuations from the physical regime. We extend our terminology 
by referring to these cases as "physical" . 

An adjacent "near-physical" case is when p,q e D r . Then W pq (n) = 
0(1), W pq (n) = 0(k' a ), where a = 2mod(iV - r,N)/N if r is odd, a = 
2 mod (r, N) /N if r is even. 3 . The W pq (n), W pq (n) are therefore no bigger than 
W pq (0), W pq (0). If this occurs on only one or two rapidity lines (or half- lines), 
while the other (background) rapidities are as in the previous paragraph, then 
we still expect the thermodynamic limit to exist. 

The other "near-physical" case is when p g T> r , q e T> r+ \. This has the 
same behaviour, except with W pq {n) interchanged with W pq (n). Again we 
expect the thermodynamic limit to exist if this occurs on only one or two 
rapidity lines. 

We focus attention on the low-temperature case when k' is small. However, 
we do not expect any discontinuities or non-analyticities for < k' < 1, so 
expect the above remarks to generalize to < k' < 1. For instance in T>\ we 
take |/Xp| > 1. Then x p is constrained to a near-circular region enclosing the 
point x p = 1 with | arg(x p )| < tt/(2N), while y p lies anywhere in the complex 
plane except N such regions surrounding the points l,u>, . . . ,uj n ~ 1 .^^ Fig.3] 

We do need to consider whether the pictures of Figures 0J El El all corre- 
spond to the same analytically continuous function. For the physical regime 
(all Boltzmann weights positive real), if p e r , then q should be in domains 
T> r , P r _i, T> r+ i, respectively. Then the thermodynamic limits of the partition 
functions of each figure will exist. They will be analytic continuations of one 
another, but will lie on different Riemann sheets. 

Here we want to keep p, q lying in T>i for all three figures. For Figure UJ 
this is the obvious analytic continuation of the physical case. We must ask 
whether we expect the partition functions of the other two figures to converge, 
and if so whether they will analytic continuations of the first. 

The answer to both questions is yes. We start with 

p,q,heT>i , v e B ■ (45) 
3 To derive these formulae, one may have to use the alternative but equivalent form l|59|l of 18bjl 
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Consider the process outlined above for obtaining Figure 03 When one rotates 
and reverses the half-line p to the lower vertical position, it then intersects 
the background horizontal lines with edge weight functions W p ' jq ', W p i jq , where 
p' = R~ 1 p,q' = h. Thus p' e Vq, q' e T>±. This is the second of the near- 
physical cases discussed above. The thermodynamic limit will still exist. 

We then shift the background rapidities v, h one half-step down to T>q, Bq. 
The background weights remain physical, as are now those on the half-line 
R~ x p. Now those on q have weights W vq , W vq , where v e X>q, q e T>\. This is 
the same near-physical case, so the thermodynamic limit still exists. 

Finally we rotate the line R~ 1 p through one more right- angle to assume the 
horizontal position of Figure^! Then the weights on it become W vR -i p , W v y R-^ p , 
and v,R~ 1 p e T>q. This is the first of the near-physical cases, so again the 
thermodynamic limit exists. 

Thus we can go continuously from the F pq (a) of Figure 0] to that of Figure 
|SJ keeping the partition function convergent. We therefore expect both figures 
to give the same analytic function. 

A similar argument applies to the rotations and reversals necessary to go 
from Figure to Figure El Now we shift v, h one half-step up to T>\, B\ after 
reversing p to Rp. Again, the background weights remain physical, while 
those on the two half-lines are either physical or near-physical. We expect 
both figures to give the same analytic function. 

The Fourier transform ratio G pq (r) 



As in |20j . we define 

N-l 

F pq (r) = 5>™F p9 (a) , (46) 

a=0 

G pg (r) = F pq (r)/F pq (r-l) . (47) 
Then the above equations (|4*T|) - (|4*4"|) yield 

Gn P ,Rq{r) = l/G pq (N-r + l) , (48) 

G R s P ,R Sq (r) = l/G qp (N-r + l) , (49) 
G m {r) = G Rq ^ p {r) , (50) 



and from l|47|) . 



Also, 



N 

II G p*(r) = 1 • ( 52 ) 



r=l 



G Mp ,g( r ) = G p! M-i q (r) = G pq (r + 1) , (53) 
where M is the rapidity operator such that 

{a,M P ,bMp,CMp,dMp} = {a p , w _1 c p , d p } . (54) 
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A significant point that we shall use is that if q = p, then the break in the 
rapidity line in Figure 0] disappears: the two half-lines become one, and then 
they can be removed to infinity. In the thermodynamic limit the probability 
Fp P (a) is therefore independent of p, so it is also true that 

G pp {r) = independent of p . (55) 

Since F pp (a) is the probability that the central spin has value a, 

M r = (u ra ) = Fpp(a)/F p p(0) = G pp (l) • • • Gp P (r) , (56) 

so a knowledge of G pq {r) is certainly sufficient to determine the order param- 
eter M T - 

Equations (|48[) - (|53|) can be regarded as functional relations for the func- 
tions G pq (r). In [2U] we showed how they can be solved (using an analytic- 
ity assumption) for the Ising case N = 2. We also showed how an obvious 
generalization of this result to N > 2 satisfies the functional relations for 
G pq (r)G Rq,R p (r) , but is wrong. (It disagrees with the known series expansions 
for M r -) This is a salutary lesson that these relations do not by themselves 
determine G pq (r): one has to input the correct analyticity properties. 



5 Low-temperature limit 

Now we briefly consider the low-temperature limit, when k' is small, using the 
picture of Figure |1] 

As in |2()l I26j , we focus on the case when p, q e T>i , when 

x p ,x q ~l , /j ip ,i_i q = 0{k'~ 1/N ) (57) 

so y p ,y q are arbitrary, of order unity. We also take v e Bq, h e T>\. This 
includes the physical regime where all the Boltzmann weights are real and 
positive. 

The dominant contribution to Z{a) then comes from the configuration 
where all spins other than a are zero. Then 

F pq (a) = ^(JV-a)F^(JV-a)F t/l (a)^(a) , (58) 

allowing the two vertical rapidities to have different values v, v'. Since W pq (N) - 
1, we can divide the relation (|8b"|) by W pq (N) to obtain 

N _ j 

W pq (n) = M — n :JP^ • (59) 

j=n+l WX P ~ U X 1 

Using this, (|%a|) and (with k = 1), we find that F pq (0) = 1, while for 

1 < a< N- 1, 

F ™ {a) ~ mi - - M i - wj-Hp ■ m 



i=i 



This satisfies the relation (|%2*)l . 
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As expected, the background rapidities v,v', h cancel out of these expres- 
sions. Also, when q = p, F pq (a) is independent of p. From l|46jl. 

{(jn) = = 1 _ k ,2 r{N _ r)/(2iY2) + 0(A .,4 ) (61) 

in agreement with © and 1)11(1 . (The second-order terms in F p q(0) cancel out 
of this calculation.) 

For N = 3 the author has obtained series expansions for F pq (a) to order 
/c /8 .[Hl eq. 48]. 



6 A solvable special case 

We have not determined G pq {r) for general values of the two half- line rapidities 
p and q. What we have done is obtain it for the case when {a q ,b q ,c q ,d q } = 
{a p ,u}bp,c p ,d p }, i.e. 

Xq = x p , y q = ujy p , n q = fx p . (62) 

The calculation and the result are very similar to the calculation of the free 
energy of the T2{t q ) model. 

We have given the expressions (J5J) for the Boltzmann weights in terms of 
the original rapidity variables a p ,b p ,Cp, d p because they make it clear that the 
weights remain finite and non-zero when b p , b q become zero. This corresponds 
to taking y p = y q = 0. Hence the special case (|fi2*|) then intersects the desired 
case q = p. From (|55|) . this is sufficient to determine G pp (r) and hence the 
order parameter M. r . 

In this section we shall consider the three pictures in Figures 0J El El of 
the function F pq (a). They all give the same function F pq {a) in the limit of an 
infinite lattice. However, they do differ for a finite lattice. We shall use this 
difference to manifest different properties of F pq (a). 

We shall then put these properties together to make what we believe to 
be a plausible and correct assumption as to the analyticity properties of the 
Fourier transform ratio G pq (r) for < k' < 1. This assumption is the key to 
derivation of this paper. 

One function that we shall use is 

e(r) = l-NS r>0 , (63) 

where <5 rj o is to be interpreted modulo N. Thus e(0) = e(N) = 1 — N, and 
e(r) = 1 for r = 1, . . . , N — 1. 



First picture 

Consider the definition (|34fl of F pq (a), where Z{a) is the partition function 
of the lattice shown in Figure 0J For the moment ignore the restrictions (J62j) 
and allow p, q to be general rapidities. For a finite lattice (with boundary 
spins fixed to zero), it is readily seen that 

Fpq(a) = ( / u p / / u (? ) a x (rational function of x p ,y p ,x q ,y q ) . (64) 
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Consider p p ,x p ,y p as functions of the complex variable t p . From (JJJ) they 
are multi-valued functions, with N branch cuts Co, . . . ,Cj\/_l- The cut Ci is 
along the straight-line segment (uj 1 p,oj % / p), where p = [(1 — k')/(l + k')] 1 ^ . 
For iV = 3, these cuts are shown in Figure [7| On them \p p \ = 1. Here we 
concentrate our attention on the case where p eT>i, which is when 

\H P \ > 1 , (argxpl <vr/(2iV) . (65) 

For k' small, x p is restricted to a small region 1Z round x p = 1. We shall say 
x p £ 1Z, or simply x p ~ 1. On the other hand, y p can lie almost anywhere 
in the complex plane, being excluded only from 1Z and corresponding small 
regions round y p = to, . . . , u> N ~ % . 

For a finite lattice, it follows from (j64| that F pq {a) is also a multi- valued 
function of t p . To make it single-valued, we must restrict t p to the cut plane 
Figure 

Now consider the low-temperature limit of section 5. Making the substi- 
tutions (JMJ) into (Hp), we observe that there are significant cancellations. We 
obtain 

*«M = S57i 1 ^ + 0(tf 4 ) . (66) 



From (gnj), (S3) it follows that 



G(r)-1 (N ~ 1 = 2r)fc/2 £ ( r ) fe/2 | Q(t 

for r = 0, ...,N- 1. 

To this order there is no evidence of any zero or singularity in G pq (r) for 
any value of t p other than one. This singularity is consistent with there being 
a branch cut on the real axis, as in Figure but at this very preliminary 
stage there is no evidence that the other N — 1 branch cuts, or any other 
singularities, occur. 

Second picture 

Now we consider the picture of Figure |S] for Z{a) and write 

p' = R- l p (67) 

for the rapidity of the lower half-row. The rapidity of the upper half-row is q. 
Then 

x p — Up' i Up — wXp' , p p — 1/ p p i i so from (|62|) the rapidities of the 
two rows are related by 

x q = y P ' , y q = u 2 x p ' , p q = i/p p > . (68) 

For the moment, take the lattice to be finite, of half-width s. In Figure 
we have shown the sites of C immediately above and below the two half- 
rapidity lines (including the left-hand site with spin a) as open circles. All 
other sites are shown as filled circles. 

If we duplicate the spin a as two open circles at the same heights as the 
others, these open circle spins correspond precisely to the spins in Figure [21 



„/4\ 
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Figure 7: The cut t p -plane for N — 3, showing the contour C. 



The spins between them lie at the centre of stars, as in Figure El Summing 
over them, we obtain precisely the weight (|19|) . with a\,a[ = a. Because the 
boundary spins are fixed to zero, we also have o~ s , a' s = 0. 
Comparing (|68j) with we see that we have 

k = , t = 2 , j = 2 , (69) 

so we can use the results of the j = 2 sub-section of Section 3 for the contribu- 
tion to each term in Z(a) of edges crossing the half-lines R~ l p, q. Thus Z(a) 
is a sum of expressions ()19|) (with coefficients independent of p, q), so from 
((2*3*)) Z(a) is a polynomial in t p of degree s — 1. It has no other dependence 
on p. 

Again we consider the low-temperature case (|57[) . The dominant contribu- 
tion to any Z(a) comes from the configuration where all the remaining spins 
other than a are zero. Also, the other four edges incident to spin a give a 
weight that is unity when a = 0, while when a ^ to leading order their 
contribution is k' 2 /[N 2 (l - uj a )(l - uj~ a )}. 

Using Table [*Q with p replaced by p' and y v = y v i = 1, and noting from 
the equations of this sub-section that ut p i = t p , it follows that 

z(o) = (i-tpY^ + oik' 2 ) , 
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and for a ^ 



7i \ k> 2 (i-uj a t p )(i-t p y- 2 4 

Again we obtain the leading-order result 1)66(1 for F pq (a), but now we see also 
that, for s finite and 

F w (a) = Z(a)/Z(0) , (71) 

where the -Z(a) are polynomials in t p of degree s — 1, equal to (1 — when 
fc' = 0. 

Hence, by continuity, for sufficiently small k 1 the zeros of the polynomials 
Z(a) must lie close to t p = 1. 

As we take the limit of s large, we only expect singularities to occur in the 
vicinity of these zeros. This reinforces the suggestion of the last sub-section 
that in the limit of a large lattice, F pq {a) and G pq {r) may be analytic functions 
of t p except for singularities in some region surrounding t p = 1. Again, there 
is no evidence of the branch cuts in Figure other than the cut on the real 
positive axis. 



Third picture 

Finally we consider the picture of Figure El for Z(a), and write 

P =Q , q =Rp (72) 

for the rapidities of the half-rows immediately below and above the spin a. 

Then x q > = y p ,y q > = u>x p ,fi q / = l/fJ, p , so from (|62j) the rapidities of the 
two rows are related by 

Xq i=u)~ l y p i , y q = Lox p < , n q = l/n p > . (73) 

Again we initially take the lattice to be finite, of half-width s. The open 
circle spins in Figure |HJ other than a, correspond to the spins in Figure |3J 
Let the lower such spins be a±, . . . , a s , the upper ones a[, . . . ,a' s . Summing 
over the spins between them, we obtain the weight (|19|) . with the boundary 
spins cr s ,a' s = 0. We also obtain the weight of the two edges incident to a. 
From and l|17j). these contribute a further weight V vq f p '(— a[, —o\\ —r) to 
F pq (r). Using (|24|). we see that the total factor contributed by the edges that 
cross the lines q, Rp is J r (ci, . . . a' s ). 

From l|75|>. 

k = -1 , l = N + 1 , j = N , (74) 

so we can use the results of the last sub-section of section 3, with v = N — 1. 
Summing over all the spins other than a in Figure El we see from (|25jl that 

Xp {r) F pq (r) = Ar/Ao , (75) 

where 

m(r) = mod(iV -r- 1,N) (76) 
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and each A r is a weighted sum of functions J T , with coefficients independent 
of p. Thus A r is a polynomial in t p of degree (N — — 1). This is the 
complete dependence on p. 

When k! = 0, then x v = x v i = 1 and the only contribution to the sum- 
over-states is from all spins being zero, so from (jHSJ) 

JV-1 

A r = IKl-cvHpY- 1 . (77) 

From continuity, for sufficiently small k' we see that all the zeros of the poly- 
nomials A r must be clustered round the iVth roots of unity in the complex 
ip-plane, excluding t p = 1. 

So now we have information that is reciprocal to that of the previous sub- 
section. In the thermodynamic limit of s large we expect x™^ F pq (r) to have 
singularities near the t p = to, to 2 , . . . , u;^ -1 , but not near t p = 1. From (J75|) 
and (jUJ), 

G pq (r) = xf^ArlAr-i , (78) 

for r = 0, . . . , N - 1. 
The relation (I5T1). 

So far we have not used any of the six functional relations (|48p - ()53[) , Half 
of them are not helpful for the present argument, since if the arguments p, q 
of the function G on one side of the relation satisfy our restriction 1)621) , then 
those on the other side do not. 

One relation that is helpful is ()51|) . Set 

Pl = R-\ , Ql = Rp , (79) 

then, using (|35|) and (|62|) . we obtain 

so pi and q\ also satisfy the restriction l)62() . The relation ()51|) becomes 

G » M = ^5i G -« M - < 81 > 

which implies 

^Gpqir) = V p X G Pim {r) (82) 
for r = 1, . . . , N — 1. Also, from ()52l) it follows that 

Xp e{r) G pq (r) = y p e{r) G Puqi (r) (83) 

for r = 0, . . . , N - 1. 

If we write p = {x p , y p , fi p }, then p\ = {y p , x p , 1//U P }. Hence p\ is obtained 
from p by interchanging x p with ?/ p and inverting which leaves t p inter- 
changed. This is what one obtains by taking t p across the branch cut on the 
positive real axis in Figure [7| and then returning it to its original value. The 
relations (|82|) , (|83|) therefore say that the functions on the lhs of the equations 
also return to their original values, i.e. that they do not have this branch cut 
in the complex t p plane. This agrees with (|78|) and our observations that we 
do not expect the ratios of the A r to have any singularities in the vicinity of 
t P = 1. 
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7 Analyticity assumptions 



Define, for r = 1, . . . , N, 




= x, 



p 



t{r) G pq (r) 



(84) 



The functions G pq {r) are defined by (jjEj) and (|4T|) in terms of F pq (a). 
This in turn is proportional to Z(a), the partition function of the lattice with 
the broken rapidity line. Three equivalent forms of this lattice are shown in 
Figures El El El 

In the previous section we took q to be related to p by (|d^)) . We presented 
evidence based on Figures 0] and El to suggest that for sufficiently small k' the 
functions G pq (r) are analytic and non-zero in the complex t p -plane, excepting 
some region surrounding and not far from the point t p = 1. Then we used 
Figure El to obtain evidence suggesting that the functions H pq {r) are analytic 
and non-zero except near t p = u>, . . . , uj n ~ 1 . We showed that this is consistent 
with the functional relation ()51l) . 

We remarked at the beginning of section 5 that the chiral Potts Boltzmann 
weight functions W pq (n), W pq (n) remain finite when b q = and a q ,c q ,d q ^ 0. 
For this reason we neither see nor expect any singularities in the G pq (r) at 
t p = 0. The same is true when c q = 0, implying that the G pq (r) remain finite 
and analytic when t p — > oo. 

We therefore make the following analyticity assumptions, taking p, q to 
be related by (|62|). and regarding G pq {r),H pq (r) as functions of the complex 
variable t p . We take both p and q to lie in the domain Pi, i.e. to satisfy 
(|65|). Although our evidence is obtained when k' is sufficiently small, we do 
not expect any non-analyticities in k! throughout the ferromagnetic regime 
< k' < 1. We therefore expect the assumptions to be true for < k/ < 1. 

ASSUMPTIONS 

There exists a closed contour C surrounding the branch cut on the positive 
real axis in Figure such that 

i) the functions G pq (r) are analytic, bounded and non-zero outside and on 



ii) the functions H pq (r) are analytic and non-zero inside and on C. 

The functional relation ()51|) implies that H pq (r) has no branch cut on the 
positive real axis, so (ii) follows if we strengthen (i) to apply to the whole 
tp-plane except for the cut on the positive real axis. 

These assumptions are very similar to those we made to calculate the free 
energy of the T2(t q ) model by the "inversion relation" method |221 p. 419]. In 
fact we shall find that G pq (r) has a similar form to that free energy. 

These analyticity properties can be expressed as symmetries. |221 12"7] For 
instance, if we exhibit the dependence of G pq (r) on x p ,y p by writing it as 
G r (x p ,y p ), then Assumption (i) implies 



C. 



for J = 1 
Pi. 



G r (x p ,y p ) = G T (u 3 y v ,Lo 3 x p ) , 
. . . , N — 1, (x p , y p ) g Pi and (ujiy p , io~'- l x p ) in a domain adjacent to 
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8 Calculation of G pq (r) 

The above assumptions are sufficient to calculate each of G pq (0), . . . , G pq (N — 
1) to within some constant factor. These factors can then be determined from 
the property l|55|). 

The we can write 1)84(1 as 

x f) = Gpq {r)/H pq {r) . (85) 

Then it follows immediately from the above Assumption that we can obtain 
the functions G pq (r), H pq (r) by a Wiener-Hopf factorization of x p . (See for 
example eqns. (50) - (53) of [H].) 

The variables x p , [i p are function of t p : let us write them as x(t p ), n(t p ). 
The function log x(t p ) is single- valued and analytic on C in the i p -plane. Define 

where the integration is now round the corresponding curve C in the complex 
s-plane, and t lies outside C. Define B + (t) by the same equation, but taking 
t to be inside C. In the first case shrink C to a curve C_ just inside C, in 
the second expand it to C + lying just outside C. Now take t to lie between 
C+ and C-. Then the combined curve of integration for B+(t) — B-(t) can 
be deformed to a closed curve surrounding t and close to C, so by Cauchy' 
integral formula 

fl+(t)-S_(t) = logx(t) . (87) 

Hence from (|85[). 

e(r)5_(t) + logG M (r) = e(r)B+(t) + log H pq {r) . (88) 

From (|86|) and our assumption, the lhs is analytic and bounded outside 
C. Similarly, the rhs is analytic and bounded inside C. Each side is therefore 
entire and bounded, so be Liouville's theorem each is a constant. Hence 

logG pq (r) = constant — e(r)B^(t p ) , (89) 

where the constant is independent of p, but may depend on r. 

This result is of course consistent with our assumption, and with the 
stronger assumption that G pq (r) is analytic and non-zero in the t p -plane, ex- 
cept for the cut along the real axis in Figure [7| The contour C can be shrunk 
to just surround this cut. 

Integrating by parts, we can write (fS^|) as 

B_(t„) = = --^ / log(s - t p ) ^4 ds > ( 9 °) 
2-7T1 Jc x(s) 

The quantity ^(s) N moves once around the unit circle as s moves around 
C (in the positive directions). Setting n(s) = exp(—i9/N) and using (j%9*)) . we 
obtain for r = 1, . . . , N, 

G pq (r) = C r S(t p y^ , (91) 
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where C r is a constant, 

log S(t p ) = - A log k + J_ £ J^- d iog[A((9) - t p ] d8 , (92) 

and 

A(0) = [(l-2k'cose + k' 2 )/k 2 ] 1/N . (93) 

The log A; term has been included in (|92j) to ensure that S(0) = 1, choosing 
the logarithms to ensure that S(t p ) is positive real when t p is real and off the 
branch cut. The condition (|52|) is satisfied provided that 

dC 2 ---C N = 1 . (94) 
Calculation of Ci, . . . , CV- 

We still need to determine the constants C\,... , Cjv- We noted at the be- 
ginning of section 5 that when y q = y p = 0, with the restriction (|62[). 
G pq (r) = Gp P (r). In this case t p = and we can evaluate the integral in 
(EH) , giving 5(0) = 1 and 

G pp (r) = C r . (95) 

Now consider the case when a q , . . . ,d p satisfy the restrictions given imme- 
diately before ()62j) . and c q = c p = 0, so y q = y p = cxd. It is not true that 
q = p, but from (|54|) it is true that q = M~ l p. Now t q = t p = oo and the 
integral in (|92j) is zero. Using H53|) . we therefore have 

G pp (r + 1) = G PiM -i p (r) = C r k~ 2 ^' N2 . (96) 

Because of (|53|) . the G pp (r) in these last two equations are the same. 
Eliminating them, we obtain 

C r _x = /t 2 ^ 2 C r (97) 
for r = 2, . . . , N. From these equations and it follows that 

C r = A;^ 1 - 2 ^ 2 (gg) 

for r = 1, . . . , N. 

This completes the calculation of G pq (r) for the case when p, q are related 
by (jH2j): it is given by (|91|) and (|98|) . 

We also immediately obtain from (|95|) that 

G pp (r) = fc( Ar+1 - 2r )/ iv2 (99) 

for r = 1, . . . , N . Hence from (|5fij). 

A^ r = = k r{ - N ~ r ^ N2 , (100) 

for r = 0, . . . , N. This verifies the conjecture ©. 
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Relation to T 2 (t p ). 

The free energy of the T 2 (t p ) model is given in eqn. (39) of jl()j and (73) of 
|22| . These expressions are very similar to the integral in 1)91(1 . In fact if we 
use eqn. (73) of [22] for T2(t p ), and manifest its dependence on X q by writing 
it as T 2 (t p , X q ), then 

k l l N S{ut p f N = r 2 (t p ,0)/r 2 (t p ,k') . (101) 
The product G(t p , r) ■ ■ ■ G(uo N -%, r). 

Exhibit the dependence of G pq (r) on t p by writing it as G r (t p ). Consider the 
product 

L r (t p ) = G r (t p )G r {ujt p ) • • • G r {u N -H p ) . (102) 
This can be evaluated from (|91|) . giving: 

L r (t p ) = D r x e } r) , (103) 

where 

D r = C?k~< r ^ N = k l+& ^- 2r ' N (104) 

for r = 1, . . . , N. 

In fact we could have obtained this result very directly from the assump- 
tions of section 6. These imply that the functions x p e ^G r (t p ), G r {ujt p ), ... 
,G r (u N ~ l t p ) are all analytic on and near the positive real axis. Hence so is 
the product L r (t p )/x e p r \ Since this is unchanged by multiplying t p by u, 
none of the branch cuts in Figure El appear, so it is analytic everywhere. It is 
bounded, so by Liouville's theorem it is a constant. This verifies (|103[1 . but 
does not determine the constant. To do this, repeat the argument of (|95() - 
(|9Hjl. From 0, x p = k l l N when y p = 0, and x p = k~ l l N when y p = oo, so we 
obtain 

G^rf = D r k< r ^ N = D r ^k-< r - l ^ N . (105) 

From DiD 2 ■ ■ ■ D N = 1. It follows that D r is indeed given by (fTUIj) . 

It also follows that G pp {r) is given by (|99]> . which verifies the conjecture (j3J). 
This route avoids using the Wiener-Hopf factorization and the integral ((911) . 



9 Other special cases. 

It is natural to ask if we can evaluate G pq {r) when p, q are related, not by 
1)62(1 . but by the more general condition 

x q = x p , y q = to l y p , n q = \i v , (106) 

where i = 0, . . . , N — 1. 

For such a case, let us introduce an extra index i and write 

jV-l 

G pq (r) = G ir {t p ) , L ir (t p ) = J] G ir (Jt p ) . (107) 

3=0 
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As a first step, one can use the series (48) of to expand L ir (t p ) to order 
k' 8 for N = 3. Indeed, this is how the author discovered that the case (|62|) 
might be tractable. For i = 2 and r = 0, 2 we have not observed anything 
particularly simple, but for i = 2, r = 1 the series are consistent with the 
conjectures 

L 21 (t p ) = x 2 p , G 2 i(t„) = k 2 / 9 S(t p )S(iot p ) , (108) 

<S(t p ) being the function defined by (|92j). 

We can use the general result (|qT?)) to expand Gi r (t p ) to first order in k' 2 . 
For i = 0, . . . , iV — 1 and r = 1, . . . , N — i the results simplify, giving 

0"'(« = i + ^(* + *-"-i-»tOT5) • 

which is consistent with the formula 

G ir (t p ) = k { - N+l ~ 2r ^ N2 S{t p )S{ujt p )---S{uj i ~ 1 t p ) . (109) 

This formula also agrees with (|99|) when t p = 0,q = p, and when i p = oo, q = 
M~ % p. From preliminary calculations we suspect (|109fl can be justified in 
the same way that we justified Assumptions (i) and (ii), by arguing that 
G pq (r) has no branch cuts on C\, . . . , Cjy-i, while H pq (r) has no branch cuts on 
Co,C_i, . . . ,C\-i. We could then derive ()109|) by a Wiener-Hopf factorization, 
taking the contour C to surround the cuts Cq,C_i, . . . ,Ci_$. 

For other values of i,r awkward factors such as 1/(1 — cj) occur that 
multiply the 1/(1 — up tp) terms, and the conditions to use (|21|) in (|24| are 
no longer satisfied because n = —r is not in We have observed no simple 
patterns for these cases. 

10 Summary. 

We have derived the long-conjectured result (jSJ), (|11|) for the order parameters 
of the ferromagnetic chiral Potts model. The method does depend on the 
analyticity assumptions of section 7, but in this respect it is no different 
from the standard derivation of the chiral Potts model free energy. In both 
calculations one argues that a certain u T2(t p )" type-function is analytic except 
for a single branch cut, and obtains a rule relating the two values of the 
function on either side of that cut.p H? H 122] 

In one respect this order parameter calculation is easier than that of the 
free energy. Although we have performed a Wiener-Hopf factorization to 
obtain the function S(t p ) of section 8, we did not need to. As we remarked 
at the end of that section, it is sufficient to calculate the elementary function 
L r (t p ). 

It seems that it is this function that is the desired generalization to arbi- 
trary N of the simple N = 2 function L pq {r) of [20] • For N = 2 the chiral Potts 
model reduces to the Ising model. There is a difference property and all func- 
tions depend on p, q only via k and u q — u p , where lip, Uq dbTQ. elliptic function 
arguments. For a long time the author thought one had to keep p, q indepen- 
dent to avoid the trivial situation when u q — u p is merely a constant. Of course 
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this is not so: for N = 2 the restriction (|62|) corresponds to u q = 2K — u p , K 
being a constant (the elliptic integral). Hence u q — u p = 2K — 2u p and the 
functions do not merely degenerate to constants. It is this " superintegrable" 
N = 2 case that we have generalized in this paper. 

The actual calculation is not difficult, being given in section 8. Most of this 
paper is concerned with presenting arguments for the analyticity assumptions 
of section 7. 

For N = 3, low-temperature expansions of G pq (r) were developed in 
As we remarked above, these were very useful in developing and checking the 
ideas that led to our analyticity assumptions. They agree with the results 
iPTTjl. (PJi. and with the conjecture (|T09|) . 

We have only obtained the generalized order parameter function G pq (r) 
when p, q are related by (|1U6|) . with i = and 1. At the end of the previous sec- 
tion we have also conjectured how these results may extend to i = 2, . . . , N — 1, 
r = 1, . . . , N — i. It would be of some interest to study the other cases of (|1U6|) . 
and indeed the general case of no restriction on p, q. This last is not an easy 
problem: we intend to comment on the difficulties in a subsequent paper. 
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